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We derive a detection method for a stochastic background of gravitational waves 
produced by events where the ratio of the average time between events to the 
average duration of an event is large. Such a signal would sound something like 
popcorn popping. Our derivation is based on the somewhat unrealistic assumption 
that the duration of an event is smaller than the detector time resolution. 

Consider a large collection of similar gravitational wave sources. If we cannot 
resolve the individual signals produced by these sources and know only their sta- 
tistical properties, the signals form a stochastic background. We group stochastic 
backgrounds into the following two classes: those for which the individual sources 
have (i) average redshift z 3> 1 and (ii) average redshift z ~ 1. Typically in case 
(i) the number of individual sources is large enough that the duration of events 
is long compared to the time between events, and then the central limit theorem 
enforces Gaussianity of the background. However, for some sources in case (ii) the 
anticipated number of individual sources might be small enough that this argu- 
ment breaks down and a non-Gaussian stochastic background might be expected. 
Examples of (ii) are supernovaeu, binary inspirals, and young neutron stars. Well 
developed data analysis methods exist for detecting Gaussian backgrounds. Here 
we derive a data analysis method for a non-Gaussian background under idealized 
assumptions; further work is needed to obtain realistic detection methods. 

Assume a pair of co-located, aligned detectors with outputs h\ = n\ + s k (i = 
1,2) of uncorrelated Gaussian noise n k and possibly a stochastic background signal 
s k . Here i labels the detectors and k with 1 < k < N labels the time samples 
of the detector output. We assume the following probability distribution for the 
stochastic background signal: 
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where s = (s 1 , . . . , s ). Thus, the individual time samples are assumed to be 
statistically independent and identically distributed; this is because we assume 
the duration of events is smaller than the detector sampling time, so each event 
affects only one time sample^. In Eq. (Q), the parameter £ is the probability that 
a randomly-chosen data sample will contain an event, and a is the signal strength. 
The Gaussian case is recovered in the limit £ — > 1. 

We use Bayesian methods to motivate the choice of detection statistic. Let 
P' l > be the probability that a signal is present in the data (posterior) and P^ 
be the probability that a signal will be present (prior) . The relation between these 



a Detector time resolutions can be on the order of milliseconds and even supernovae waveforms 
can be on the order of hundreds of milliseconds, so this assumption is somewhat unrealistic. 
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probabilities is given in terms of the likelihood ratio A by 

p(l) p(0) 



l-PW A 1-PW ^ 

A natural detection criterion is to threshold on which is equivalent to thresh- 
olding on A. Thus, we will say that the signal has been detected if A > A* for some 
threshold A*|]. 

To compute A we need only detector data and the probability density functions 
for the signal and for the noise in each detector. We find 

A = / d N s PM^(hi- s,h 2 -s) (o) 

where is the vector (h\, . . . , hf), pA/"iJV 2 ( n i? n 2) is the joint probability density of 
the detector noises ni and n 2 , and Pg\s) is the probability density that the signal 
s is present in the data, given that some signal is in the data. We take pg (s) = 
/ df J da p(s|f,a)p(°>(f,a), where p(°>(f, a) is the prior distribution for the signal 
parameters £ and a. Combining these equations and assuming uncorrelated white 
Gaussian noise of unit variance in each detector^] we obtain 



A = J d£ J daA(^a)pW(^a) (4) 

where 

baA(£,a)=^ln(l + £0 (5) 
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The detection statistic we suggest is the maximum likelihood statistic 

A max = max A(£, a), (7) 

a, 5 

which we expect to be large when A itself is large. The maximization can be 
carried out numerically^. If a signal is detected, then the values £ max and a max 
which achieve the maximum give estimators of the true values of a and £. One can 
also form a posterior distribution for a and £ from the function A(a, £) and thereby 
attempt to distinguish between Gaussian and non-Gaussian backgrounds. 
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k We restrict attention here to finding a useful statistic and ignore the issue of how to choose the 
value of the threshold, which may be done using either Bayesian or frequentist methods. 
c We expect to be able to relax this assumption in the near future. 

d While such a search is computationally intensive, we have had some success with numerical 
simulations using small (length ~ 10 6 ) data segments. 
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